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The energy levels of two and three anyons in a two- 
dimensional parabolic quantum dot and a perpendicular mag- 
netic field are computed as power series in 1/|J|, where J is 
the angular momentum. The particles interact repulsively 
through a coulombic (1/r) potential. In the two-anyon prob- 
lem, the reached accuracy is better than one part in 10^ . For 
three anyons, we study the combined efi'ects of anyon statis- 
tics and coulomb repulsion in the "linear" anyonic states. 

PACS numbers: 03.65. Ge, 74.20. Kk 



I. INTRODUCTION 

Recent developments in semiconductor technology 
(e.g. MBE and electron lithography) opened the pos- 
sibility to create totally confined electron. systems, the so 
called artificial atoms or quantum dotaj'D. This is one 
of the examples of present-day-physics' interest in con- 
fined finite systems, among whicih one can mention also 
the atomic and electronic trapgj, and the condensation 
of confined bosonsQ. 

Quantum dots exhibit very interesting properties like 
the possibility of varying their parameters (number 
of electrons, applied fields, dot's geometry, tempera- 
ture) over a. wide range, the observation of conductance 
oscillationsO and Kohn's theoremO, etc. 

Theoretically, these systems have been studied mainly 
with the help of numerical methods. However, analytic 
approaches have proven to work extremely well provid- 
ing, at the same time, a qualitative understanding of the 
quantum dynamics. Among these aoproaches one can 
mention the semickissical quantisationQ, regularisec 
turbation theoryolj, Padej-approximant techniquest 
and the 1/A^-expansionu llJ. 

In the present paper, we continue the analytic- 
qualitative line of research and apply the 1/A^-expansion 
{N is the absolute value of the angular momentum) to 
compute the energy levels of two and three anyons in a 
model parabolic dot. The particles interact through a 
coulombic {1/r) repulsive potential. A magnetic field is 
applied perpendicularly to the plane of motion. 

Numerical results for the two-anyon system were ob- 
tained by Myrheim et alllEl. Exact analytic solutions at 
pacticular values of the coupling constants were found 
irJlB. Bohr-Sommerfeld quantisation was applied to this 
system at low magnetic fieldsll3. We shall see that our 
method provides extremely accurate solutions for states 



with angular momentum \J\ > 2. A picture for the "ge- 
ometry" of the states (the spatial distribution of prob- 
ability) is obtained also. In the three-anyon problem, 
however, to our knowledge there are no numerical or ap- 
proximate calculations. 

We start from the Hamiltonian of Na anyons moving 
in a two dimensional quantum dot in the presence of a 
perpendicular homogenous magnetij£_field. In the bosonic 
gauge, it is given by the expression!!^ 
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in which the vector potential is taken in the symmetric 
gauge, 



Aj = -B X fi, 
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(Ti is the statistical vector potential, 
n X (rj- - fi) 
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n is the unit vector perpendicular to the plane of motion 
of the anyons, e is the anyon's charge, v is the anyonic 
parameter, wq is the frequency of the parabolic potential 
needed to confine the anyons in the dot and e is the di- 
electric constant of the medium. A dimensionless Hamil- 
tonian is obtained by means of the change of variables 
fi -^ y/h/{mCl) n 
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where Uc — eB/(rnc) is the cyclotronic frequency, fl^ = 
uj'^/A + ujq, and /S'^ = \/me^/{e^rth^) is the square root 



of the ratio between the coulombic and oscillator charac- 
teristic energies. The problem has one exactly solvable 
limit: a low-density limit, which we call the Wigner limit, 
reached when /3 ^ cxd. In the /3 ^ (oscillator) limit, 
the two-anyon problem is exactly solvableEj, whereas the 
three-jaHiyon system has an infinite family of exact linear 
statef£j. In real semiconductors, /3 ^ 1. 
Introducing Jacobi coordinates. 



Pk 



PNa 



2k 



1 






rk+1 



l<k<Na-l, 



1 






^S' 



(5) 
(6) 



the centre-of-mass and relative motions are separated 
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where pi — \pi\. We have "renormalised" the coupling 

constant, (5'^ — (3'^ / \J\^ , and the statistical parameter, 
v = V I J ^ in order to take account of the Coulomb repul- 
sion and the statistical interaction in a nonperturbative 
way when taking the formal limit |J| — > oo. We shall 
look for symmetric eigenfunctions of /i, i.e. \J\ shall be 
even. 

In the I J| — > oo limit, the only term surviving in the 
Hamiltonian is the effective (classical) potential energy 
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is the centre of mass Hamiltonian and 
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is the Hamiltonian of the relative motion. Wc introduced 
the following notation: r^j = fi — f'j, and plj = pi — pj. 

We will obtain approximate expressions for the energy 
eigenvalues of Hrei/i^^) for Na = 2 and Na = 3 as a 
function of /3 by means of the 1/| J|-expansion. 



II. THE TWO-ANYON SYSTEM 

In the two-anyon problem, we have only one Jacobi co- 
ordinate, pi, and the Hamiltonian of the internal motion 
reads 

Hrel 2 , ^c ^ . ^ ^ ^ ^ , „ ^ Pi X pi 1 2 
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Notice that n ■ [pi x pi) = J. After the scaling trans- 
formation pI ^ \J\ pI, we get 



Minimising U^/f, we obtain the leading contribution 
to the energy, eg — Ueff{poi), where the radius of the 
"Bohr orbit" is obtained from 
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Substituting pi = poi + yi/\J\^^'^ in the r.h.s. of (11) 
and expanding, we get 



E' 



-|J|V2' 

where the operator coefficients are given by 
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Similar series are written for the wave function and the 
scaled energy, that is, 
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Inserting (19), (20) and the series expansion for the 
Hamiltonian into the Schodinger equation, we may com- 
pute the coefficients tpi and e; in a systematic way. Up 
to second order, for example, the system is described in 



terms of small oscillations around the equilibrium orbit, 
i.e. the wave function is 



*o = e^'' I «), 
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where 6 is the angle associated to the vector pi, the | 
n) are two-dimensional harmonic oscillator radial states 
with frequency wi = ^^3 + 4(i' + l)^/^^^, and the first 
two coefficients for the energy are 
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Afterward, we may take account of anharmonicities. 
The results for the next two coefficients are the following 
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Notice that in both Wigner (/? ^ oo) and oscillator 
(/3 -^ 0) limits the cpcrections e^ and ee go to zero. The 
expressions found inta are reproduced if we take v — Q. 

We show in Fig. 1 the relative weight of eg in e for 
the first states with 3 = 1 and 6. The parameter v 
was fixed to 1/2 (semions). The relative contribution of 
eg is never greater than 5 x 10~^ or 3 x 10~^ for J = 
2 or 6 respectively. This shows that the l/|J|-series is 
extremely well behaved. i— , 

A comparison with the exact solutions found int£l 
is carried on in Fig. 2, where the relative difference 
|e — ee^actl /cexaci is plotted against v. The state with 
J = 6, n = is shown. It may be easily verified 

that i^e^act = pi^+"'(l+Pl/v/2|./+H + l)e"''', ^e^act = 

I J + 7/| + 2, are exact solutions of the two-anyon problem 
at 0^ = •\/2| J + 7^1 + 1. The comparison shows that the 
relative error of our estimate is not greater than 10~*. 



III. THE THREE-ANYON SYSTEM 

The internal Hamiltonian of the system of three anyons 
in Jacobi coordinates p\ and pi is written as 
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Doing the same scaling transformation p\ — > | 
and making explicit the dependence on | J|, we get 
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where pi = |pi|, p2 = |p2|, and cos6' = pi • p2/(piP2)- 
The "renormalised" /S"^ = /?'^/|Jp and v = v/J were 
introduced. 

The minimum of the classical potential entering (27) 
is reached in the configuration of an equilateral triangle 
(poi = P02, d = ±7''/2). We choose, for example, poi = 
P02, = 7r/2. poi is obtained as the solution of the 
equation 
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Then, introducing pi = poi + ?yi/vkl, P2 = Poi + 
2/2/\/PT and 6* = 7r/2 + ^/i/pT in the r.h.s of (27), we 
obtain for the Hamiltonian h a series hke (14). The first 
operator coefficients are given by 
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The operator /i2 is diagonalised by changing variables 
2/1 = (2/s + 2/m)/A/2, J/2 = (2/s - 2/„i)/\/2, z = \/2z„/poi, 
and making the "gauge" transformation h'2 = e*^/i2e~*' 
where / = sign(J)2/m2;m/(2pQ]^), the results is 
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where hg describes the motion of a harmonic oscillator in 
the coordinate j/s (the symmetric mode), 
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with wi — 1/3 + {iv + 1)^//9q]^, and ft-m accounts for two 
dimensional motion in a "fictitious" magnetic field (the 
mixed mode) 
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where W2 = ^3/2 - (3i> + l)2/(2p^i). 

The first two coefficients in the expansion for the en- 
ergy are 
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e2= oJiiria + -\ + UJ2 {2n + |m| + 1) + CJ3 sign( J)to, (38) 

where lo^ — {3iy + l)/poi7 and the quantum numbers n^, 
n, and m may be used to approximately label the states. 
Up to this order, the wave function is given by 



I'd = e''^" I J,ns,n,m), 



(39) 



where S accounts for overall rotations of the system, and 
I J, Us, n, m) are the eigenfunctions of ft,2- We notice that, 
when /? — > the energy becomes 



Eo = \J + 3iy\ + 2 + 2ns + 2n+ \m\ + sign( J)m. (40) 

These are the "hnear" three-anyon states. We stress that 
they are obtained as harmonic excitations against the 
equilateral triangle configuration, and are not necessarily 
related to a cigar- like shape of the wave function (pi >> 
P2) as described inO. 

The set of numbers {J,ns,n,7n} compatible with the 
symmetry constraints (the wave function shall be sym- 
metric) are obtained upon comparison with harmonic- 
oscillatcp-jp|Ve functions at v — 0, (3 — 0. Details may be 
found inES'ta. An additional requirement is that the state 
at /3 = should be a linear state. For example, the low- 
est linear states are the following: |0, 0,0,0 > (the g.s., 
starting from £^0 = 2 at the bosonic end), |0, 1,0,0 >, 
and |2, 0,0,-1 > (starting from Eq = 4), |3, 0,0,0 > 
and |1, 0,0,1 > (starting from Eq = 5), |4, 0,0, -2 >, 
|2, 1,0,-1 >, and |0, 1,1,0 > (starting from Eq = 6), 
etc. The lowest state with J < is | — 6, 0, 0, >, which 
starts at Eq = 8. Of course, states with small values of 
I J| can not be described within our method. 

In what follows, we restrict the analysis to levels with 
quantum numbers n^ = n = m = 0. This leaves only the 
linear anyonic states with J = 3A;, where k is an integerE3. 
The geometry of the state is an equilateral triangle. It 
can be seen from (29) that the side of the triangle in- 
creases with V when J > 0, and decreases when J < 0. 
Thus, the coulomb repulsion is much more stronger for 
J < states, and the ordering of levels may dramati- 
cally change as /? is increased. On the other hand, for 
/? ^- 00 the side grows like poi ~ 3^'^P and becomes in- 
dependent of the anyonic jparameter, as one expects. A 
strong coupling expansiontd shows that the leading con- 
tribution to the energy (potential energy) is ~ (3"^, the 
next corrections (quantum fluctuations) are ^ 1, the an- 
gular momentum and the statistical parameter enter the 
second order corrections, which are ^ 1//?^. 

The first anharmonic corrections to the energy are 
given by 
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where h'^ and h'^ are obtained from hs and hi by means 
of a gauge transformation, in the same way as explained 
above for /12. 

For a state with quantum numbers | J, 0, 0, 0), we get 
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Polar coordinates have been introduced according to 
C^ = 2/m + ^"rm tana = Zm/Vm- The only nonvanishing 
matrix elements entering the sum (41) are the following 
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Collecting everything, we arrive to 
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It may be checked that the corrections go to zero in 
both the Wigner (/3 — > oo) and the oscillator (/3 -^ 0) lim- 
its. 

We show in Fig. 3 the relative weight of £4 in e for three 
semions in states with J = 3 and J = 6. The numbers 
are similar to those appearing in the two-anyon problem. 
Thus, we expect a similar accuracy to this order, i.e. one 
part in 10'^ or better. 

In Fig. 4, the levels with ,/ = ±6 are drawn. P is 
increased from 0.5 to 8. Notice that the coulomb effects 
are stronger for the state with negative J, and that the 
levels become flatter (as a function of i/) as /3 rises. 

In conclusion, the energy levels of two and three anyons 
in a model parabolic dot were computed by means of the 
1/1 J|- expansion. The qualitative picture emerging from 
the 1/1 J|- expansion is that of a rigid structure (an orbit 
in the two-anyon system, an equilateral triangle for three 
anyons) against which harmonic and anharmonic oscil- 
lations are developed. The coulomb repulsion is much 
stronger for negative- J states. Comparison with exact 
particular solutions for two anyons shows excellent agree- 
ment. 
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FIG. 1. Relative weight of es in e. Two anyons in states 
with n = and 1/ — 1/2 are studied, a) J = 2, b) J = 6. 



FIG. 2. Comparison between the 1/| J|-estimate and the 
exact solution found in [16] for two anyons with J = 6, n = 0. 
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FIG. 3. Relative weight of £4 for three anyons in states with 
= 1/2. a) J = 3, b) J = 6. 



FIG. 4. |J|e vs v for three anyons in states with J — ±6. 
a) 13 = 0.5, h) 13 = 8. 
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